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ABSTRACT 

The brightest giant flare from the soft 7-ray repeater (SGR) 1806-20 was detected on 2004 December 
27. The isotropic-equivalent energy release of this burst is at least one order of magnitude more energetic 
than those of the two other SGR giant flares. Starting from about one week after the burst, a very bright 
(~ 80 mJy), fading radio afterglow was detected. Follow-up observations revealed the multi-frequency 
light curves of the afterglow and the temporal evolution of the source size. Here we show that these 
observations can be understood in a two-component explosion model. In this model, one component is a 
relativistic collimatcd outflow responsible for the initial giant flare and the early afterglow, and another 
component is a subrelativistic wider outflow responsible for the late afterglow. We also discuss triggering 
mechanisms of these two components within the framework of the magnetar model. 

Subject headings: gamma rays: bursts — ISM: jets and outflows — stars: individual (SGR 1806-20) 



1. INTRODUCTION 

Soft 7-ray repeaters (SGRs) emit short-duration (~ 0.1 
seconds) bursts of soft 7- and hard X-rays. They are be- 
lieved to be magnetars, neutron stars with long periods of a 
few seconds and dipolc fields of ~ 10 15 Gauss (Thompson 
& Duncan 1995; Woods & Thompson 2004 for a recent 
review). Only three giant flares from SGRs have been 
detected so far, the brightest of which originated from 
SGR1806-20 on 2004 December 27 (Hurley et al. 2005; 
Palmer et al. 2005; Mazets et al. 2005). The follow-up ob- 
servations of the December 27 event revealed a bright radio 
afterglow with the size, proper motion, polarization, and 
flux at different radio frequencies all evolving with time 
(Gaensler et al. 2005; Cameron et al. 2005; Gelfand et al. 
2005; Taylor et al. 2005). An obvious break in the multi- 
frequency light curves of the afterglow occurred around 
day 9 after the burst and subsequently the flux declined 
rapidly as cx t a (where a ~ —2.7) between days 10 and 25. 
But the afterglow hardly faded during the period of days 
25-33. In fact it underwent a slower decay. The observed 
temporal evolution of the source size was 6(t) cx t 04±015 
on days 7 to 10 (Cameron et al. 2005). This was followed 
by a faster increase (Taylor et al. 2005; also see Granot 
et al. 2005). The proper motion of the afterglow on days 
7-17 was negligible but subsequently became significant 
(Taylor et al. 2005). In addition, the spectrum is gener- 
ally consistent with a single power law index (F v cx c' 5 ). 
Cameron et al. (2005) found (3 = -0.62 ± 0.02 on day 7 
and then the spectrum steepened from (3 — —0.76 ± 0.05 
on day 15 to ^ = -0.9 ±0.1. 

Two different models have been proposed to explain 
these observational results. In Wang et al. (2005), we con- 
sidered an initially relativistic energetic blast wave from 
the December 27 giant flare, following the suggestion that 
relativistic fireballs may occur in SGR bursts (Huang, Dai 
& Lu 1998; Thompson & Duncan 2001). This is, in princi- 
ple, the same mechanism established for GRB afterglows 
(Meszaros & Rees 1997; Sari, Piran & Narayan 1998), i.e. 



the emission from a relativistically expanding blast wave 
that forms when the outflow from the SGR sweeps up the 
surrounding interstellar medium. To explain the light- 
curve break around day 9, we invoked a broken power- 
law distribution of the shock-accelerated electrons. This 
model appears to explain the rapid fading of the afterglow 
but cannot account for the rebrightening starting on day 
25 and peaking on day 33. 

Realizing this difficulty, some authors (Gaensler et al. 
2005; Cameron et al. 2005; Gelfand et al. 2005; Granot 
et al. 2005) proposed a subrelativistic outflow model, in 
which the outflow from the giant flare initially expanded 
in a cavity at a constant velocity of ~ 0.3c and about 7 
days later it happened to collide with a thin shell, which 
led to its deceleration by a reverse shock while the am- 
bient matter was accelerated by a forward shock. In this 
model, the rapid decline of the afterglow flux was caused 
by adiabatic expansion of the reversely-shocked electrons 
and the rebrightening was due to initial coasting and sub- 
sequent deceleration of the outflow. This model predicts 
a size 0(t) cx t in the coasting phase, which is inconsistent 
with the observations of Cameron et al. (2005). Further- 
more, one has to assume an unobserved density bump to 
ignite the fireball right before the observations started. It 
is unclear whether this could happen and what could be 
the source of the density bump 

The idea of an initially subrelativistic blast wave from 
SGRs was suggested by Cheng & Wang (2003), who ex- 
plained the radio afterglow light curve of the 1998 August 
27 giant flare from SGR1900+14. For this event, a re- 
brightening could be seen although the data were sparse 
(Frail et al. 1999). In this Letter, we propose a two- 
component explosion model to interpret the abundant ob- 
servational data for the December 27 event. In this model, 
one relativistic collimated outflow is responsible for the ini- 
tial giant flare and the early afterglow, and another sub- 
relativistic wider outflow is responsible for the late after- 
glow. Our model seems to provide a unified picture in 
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understanding the two giant flares and their radio after- 
glows from SGRs. In §2, we discuss why two components 
are needed for the giant flares and radio afterglows, and in 
§3, we fit the observed light curves at 4.86 GHz and 8.46 
GHz and the source size as a function of time. In §4, we 
summarize our results. 

2. AN EXPLOSION WITH TWO COMPONENTS 

Within the framework of the magnetar model, Thomp- 
son & Duncan (2001) proposed a triggering mechanism for 
giant flares, in which a helical distortion of the core mag- 
netic field induces a large-scale cracking of the crust and a 
twisting deformation of the exterior magnetic field. This 
mechanism is consistent with the three observed timescales 
in the December 27 giant flare (Schwartz et al. 2005). It 
seems that the energy release for this mechanism includes 
two components. First, as the crust cracks, the exterior 
magnetic field deforms, which leads to a purely magne- 
tohydrodynamic instability. Such an instability probably 
produces reconnection events and induces magnetohydro- 
dynamic waves outside the star. The dissipated wave en- 
ergy is initially locked onto the magnetic field in a thermal 
photon-pair plasma, which, under its own huge pressure, 
expands adiabatically along some surface magnetic lines 
and eventually drives a relativistic collimated outflow be- 
cause of very low baryon loading. Second, when the core 
"toroidal" magnetic field floats up to break through the 
stellar crust, reconnection of the newborn surface magnetic 
field will drive an explosive outflow (Kluzniak & Ruder- 
man 1998; Dai & Lu 1998). Because this magnetic field 
inevitably brings some crustal matter and the reconnec- 
tion occurs at about the stellar radius, the resultant out- 
flow has high baryon contamination and thus it is subrela- 
tivistic. Alternatively, even if only one exterior magnetic- 
reconnection event associated with the giant flare happens 
within the region whose size is less than the stellar radius 
above the surface (otherwise the energy release is too small 
to power the observed giant flare, because the dipole mag- 
netic field steeply decreases with increasing radius, i.e., 
Bd oc R~ 3 ), the outward photon-pair flow launches a colli- 
mated outflow with very low baryon loading, as discussed 
above for a relativistic component. In the meantime, the 
inward emission impinges on the stellar surface, and its 
part is reflected back and re-emitted from a larger solid 
angle with much heavier baryon loading from the surface 
as a subrelativistic component. 

We now discuss why two components are needed to in- 
terpret the observational data. We consider a 7-ray energy 
release of the December 27 burst, Eq <~ 3 x 10 46 ergs, near 
the surface of a magnetar with radius of Rq ~ 10 6 cm 
in a time of to ~ 0.2 s, and thus the outflow luminosity 
is at least L <~ 10 47 ergs _1 , as implied by the observed 
hard spike (Hurley et al. 2005). An extremely large op- 
tical depth means a radiation-pair outflow with an initial 
temperature of 

T °=(l^) 1/4 ^ 21 <^ /2keV ' « 

where L0.47 — L /10 47 ergs -1 , i?o,6 = i?o/10 6 cm, and a 
is the Stephan-Boltzmann constant. Assuming the baryon 
loading rate of the outflow, M, we define a dimensionlcss 
entropy -q = L / Mc 2 . 



If rj > 1, this outflow will expand relativistically as 
the Lorcntz factor T cc R and the temperature T oc 
R^ 1 (Shcmi & Piran 1990). The temperature, Lorentz 
factor and radiation luminosity at the photosphere be- 
come T f = nrm[r?,77*], T ph = T x min^/r^) 8 / 3 , 1], 
and i p h = L x min^ry/??*) 8 / 3 , 1] respectively, where 

r?* = [ J L CTT/(47rTO p c 3 i?o)] 1/4 ^ lOOL^i?^ 74 (Meszaros 
& Rees 2000). As suggested by Nakar, Piran & Sari 
(2005) and Ioka et al. (2005), the giant flare might have 
arisen from the emission at the photosphere and/or inter- 
nal shocks if the outflow is variable. The observed average 
temperature of the hard spike, T spikc = 175 ±25 keV (Hur- 
ley et al. 2005), requires r\ = 93Lq^ 3 7 -R^g^ 16 - 

If 77 < 1, the outflow will also expand under its own 
pressure, which includes the gas pressure and radiation 
pressure, P = P g + P r = n p k B T + 4erT 4 /(3c), where 
n p — M / {AitR 2 m v c) is the proton number density at ra- 
dius R. Letting P g = P r at radius R , we define the crit- 
ical baryon loading rate, M CI = l&irampR^T^ / '(3/cs) = 
1.7 x 10 29 Lq^R^/q gs _1 . If the baryon loading rate is be- 
low M cr , the radiation pressure exceeds the gas pressure 
so that the temperature decays T oc J? -1 . In this case, the 
temperature and radiation luminosity at the photosphere 
decreases to 

T ph = T R /R ph = 220L^ 4 4 Xf 6 2 M 2 - 1 K, (2) 

and 

L ph = 4nR 2 ph aTX = 2.0 x 10 32 L ,4 7 i?o,6^2~5 2 ergs" 1 , (3) 

where the photospheric radius i? p h = <Jt M / (4irm p c) = 
1.1 x 10 13 M 2 5 cm with M 2 5 = M/lO^gs" 1 . However, 
if the baryon loading rate exceeds M cr , the gas pressure 
dominates over the radiative pressure and the tempera- 
ture decays as T oc R~ 2 . In this case, the temperature 
and radiation luminosity at the photosphere are 

T ph = T (Ro/R ph ) 2 = 2.0 x 10~ 5 L^R 3 ^ M^ 2 K, (4) 
and 

L ph - 1.4 x 10 4 L ,4 7 i?^ 6 M 2 - 5 6 ergs- 1 . (5) 

Thus, the temperature and luminosity of the emission from 
the subrelativistic outflow at the photosphere in both cases 
are much less than those of the giant flare, showing that a 
subrelativistic outflow model for explaining the giant flare 
can be ruled out. 

Therefore, we conclude that an ultrarclativistic outflow 
is required by the extremely high peak luminosity with 
millions of the Eddington value, hard spectrum and rapid 
variability of the initial spike emission of the giant flares. 
Furthermore, a relativistic jet model indeed provides a sat- 
isfactory explanation for the observed light curve of the 
initial spike of the December 27 giant flare, as shown by 
Yamazaki et al. (2005). It is possible that such an out- 
flow, after emitting the hard spike, retains some amount 
of energy and then drives a blast wave when sweeping into 
the ambient medium. 

On the other hand, an obvious bump at tdoc ~ 30 days 
in the 4.86 GHz light curve of the radio afterglow for the 
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December 27 event and a possible bump at i d cc ~ 10 days 
for the August 27 event imply that subrclativistic outflows 
with initial velocities of /3 n r,oc could begin to be deceler- 
ated by the ambient medium with density of n at t dcc . 
Thus we obtain the outflow mass and kinetic energy 



M nr = 1.2 x 10^ ri0 n_ 2 (t dcc /10days) d g, 



(6) 



and 



E nl = 5.4 x 10 44 /3„ r n- 2 ( W/10 days) 3 ergs, (7) 



where n_ 2 = n/10~ 2 cm~ 3 . This outflow could also drive 
a forward shock with negligible energy loss. 

In short, there could have been two blast waves after the 
December 27 giant flare. In the following we show that an 
initially relativistic shock is responsible for the early after- 
glow and a nonrelativistic shock for the late afterglow. 

3. FITTING THE AFTERGLOW LIGHT CURVE AND SIZE 

We consider two energetic shocks after the December 27 
giant flare: one initially ultrarelativistic blast wave with an 
opening angle of 9j, an isotropic-equivalent energy of E r , 
and the Lorentz factor To S> 1, and another subrclativistic 
forward shock with an energy of E nr and an initial velocity 
of ~ 0.2c. Both shocks expand in a uniform medium with 
density of n. The initially relativistic blast wave will en- 

]/3 1/3 

tor the non-rclativistic phase at t nr = 4.0^ 44 n_ 2 days 
and is expected to play a dominative role in the early ra- 
dio afterglow, where E r ,4 4 = _E r /10 44 ergs. As suggested 
in Wang et al. (2005), the electron energy distribution in 
this blast wave is taken to be a broken power-law form 
with pi and pi below and above the break Lorentz factor 
7b during the trans-relativistic stage, 



dNe 
d-fe 



% 



-Pi 



if 7, 
if 7, 



mm < 7e < lb 
> lb, 



(8) 



where 7 m i n is the minimum Lorentz factor. Observation- 
ally, such an electron distribution is not only required by 
the spectral evolution of the afterglow (Cameron ct al. 
2005) and the synchrotron spectrum of the Crab Neb- 
ula (Amato et al. 2000) but also suggested in fitting 
two gamma-ray burst afterglows (Li & Chevalier 2001) 
and TeV blazars (Tavecchio, Maraschi & Ghisellini 1998). 
Theoretically, the recent particle simulation of a relativis- 
tic two-stream instability by Dieckmann (2005) revealed a 
broken power-law distribution. Wang et al. (2005) have 
shown that this distribution can account for the frequency- 
dependent breaks of the light curves around day 9 and the 
steepening of the radio-band spectra with time (Cameron 
et al. 2005). 

The initially sub-relativistic component may largely 
contribute to the radio afterglow emission only at late 
times, as it is decelerated and the swept-up matter ac- 
cumulates. According to this scenario, the observed 
bump time t ~ 30 days corresponds to the decelera- 
tion time tdcc of this component, so we infer (3 nr , = 
R nr /ct dcc = 0.38(i? nr /3 x 10 16 cm)(t doc /30d)- 1 , where 
Pm,o is the initial velocity of the sub-relativistic com- 
ponent and i? nr is the blast wave radius at the bump 
time. Then we can constrain the isotropic-equivalent en- 
ergy of this component from the deceleration time, i.e., 



E nr = 1.5 x 10 44 n_ 2 (i dcc /30d) 3 (/3 M , /0.4) 5 ergs (from 
equation 7) . The electron energy distribution for this com- 
ponent is assumed to be a single power-law one with index 
p. This is deduced from the spectrum form of the late ra- 
dio afterglow of the 1998 August 27 giant flare from SCR 
1900+14. 

We first carried out numerical calculations of the dy- 
namics of each blast wave based on Huang et al. (2000) 
and obtained the temporal evolution of the source's radius. 
We next calculated light curves of the synchrotron radia- 
tion from the shocked matter at different frequencies, as- 
suming that the electron energy distribution has the same 
form as equation (8) and that e e and es are fractions of the 
total energy density that go into the electrons and mag- 
netic field respectively Combining the contributions from 
both components, we performed numerical fitting to the 
light curves at 4.86 GHz and 8.46 GHz in Figures 1 and 2 
respectively. In the fitting, the physical parameter values 
are taken:' T = 100, E T = 2.3 x 10 45 ergs, 6j = 0.129, 
Pi = 2.2, p2 = 3.5, and Rb = ib/lmin = 20 for the initially 
relativistic component, and E nr = 0.75 x 10 44 ergs, F m%0 = 
1.023, and p — 2.4 for the initially subrelativistic compo- 
nent. The other parameter values are n = 0.0363 cm -3 , 
e e = 0.34, and e# = 0.23. We take the distance to the 
source, d — 9.8 kpc (Cameron et al. 2005), which is 
slightly less than used in previous works. Our fitting gives 
the total x 2 /dof = 374/(79 - 11) = 5.5 for 11 physical 
parameters. 

The evolution of the source size with time was reported 
for the December 27 event (Gaensler et al. 2005; Cameron 
et al. 2005; Granot et al. 2005; Taylor et al. 2005). For 
the two-component model, the measured size at any time 
should be dominated by the brighter component. In Figure 
3, we present the evolution of the source sizes of the ini- 
tially relativistic and subrelativistic components, denoted 
by the dashed and dotted lines, respectively. From the 
light curves in Figure 2, we see that the flux density of 
the initially subrelativistic component begins to dominate 
over that of the initially relativistic component around day 
20, so that the measured sizes should be the sizes of the 
initially relativistic and subrelativistic components before 
and after this time respectively. The solid line linking both 
components in Figure 3 indicates the transition regime be- 
tween the two phases. 

4. CONCLUSIONS 

We have proposed a two-component explosion model, in 
which one initially relativistic collimated outflow can ac- 
count for the spectrum of the giant flare and the early-time 
broken light curves at different radio frequencies and the 
slow increase of the source size, and another subrclativistic 
wider outflow for the late-time rebrightcning of the radio 
afterglow and the faster increase of the source size. In 
addition, our model is consistent with the magnetar sce- 
nario of Thompson & Duncan (2001). In this scenario, two 
magnetic reconnection events may happen near the surface 
and in the magnetosphere respectively, which give rise to 
a relativistic collimated outflow and a subrelativistic out- 
flow because of different baryon loadings. Alternatively, 
even if only one magnetic reconnection event associated 
with the giant flare happens within a small region above 
the surface, the outward photon-pair flow launches a rel- 
ativistic collimated fireball while the inward flow strikes 
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the surface and is partially re-emitted as a subrclativistic 
baryonic outflow from a larger solid angle. Finally, the rel- 
ativistic component from a giant flare like the December 
27 event may be diagnosed using future high energy data, 
as noted by Fan, Zhang & Wei (2005). 
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Fig. 1. — Modelling the 4.86 GHz afterglow of the December 27 giant flare from SGR 1806-20 in the two-component explosion model. The 
physical parameter values are given in the text. The dashed and dotted lines represent the contributions from the initially ultrarelativistic 
and subrelativistic components, respectively, while the solid line is the sum of both contributions. The data are taken from Gelfand et al. 
(2005, open circles) and Cameron et al. (2005, solid squares). 
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Fig. 2. — Modelling the 8.46 GHz afterglow of the December 27 giant flare from SGR 1806-20 in the two-component explosion model with 
the same parameter values as in Figure 1. The data are taken from Cameron et al. (2005, solid squares), Gaensler et al. (2005, open circles) 
and Taylor et al. (2005, solid circles). 
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Fig. 3. — Modelling the size evolution of the radio afterglow of the 2004 December 27 giant flare from SGR 1806-20 in the two-component 
model with the same parameter values as Figure 1. The dashed and dotted lines represent the size evolution of the initially ultrarelativistic 
and subrelativistic components, respectively. The solid line is the size of the brighter one between both components. The data are taken from 
Cameron et al. (2005, solid squares), Granot et al. (2005, open circles), and Taylor et al. (2005, solid circles). 



